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Abstract: The Low Frequency Instrument (LFI) is an array of pseudo-correlation radiometers on 
board the PLANCK satellite, the ESA mission dedicated to precision measurements of the Cosmic 
Microwave Background. The LFI covers three bands centred at 30, 44 and 70 GHz, with a goal 
bandwidth of 20% of the central frequency. The characterization of the broadband frequency re- 
sponse of each radiometer is necessary to understand and correct for systematic effects, particularly 
those related to foreground residuals and polarization measurements. 

In this paper we present the measured band shape of all the LFI channels and discuss the meth- 
ods adopted for their estimation. The spectral characterization of each radiometer was obtained 
by combining the measured spectral response of individual units through a dedicated RF model 
of the LFI receiver scheme. As a consistency check, we also attempted end-to-end spectral mea- 
surements of the integrated radiometer chain in a cryogenic chamber. However, due to systematic 
effects in the measurement setup, only qualitative results were obtained from these tests. The mea- 
sured LFI bandpasses exhibit a moderate level of ripple, compatible with the instrument scientific 
requirements. 

KEYWORDS: Instruments for CMB observations;Space instrumentation;Microwave 
radiometers; Spectral response. 
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1. Introduction 

The objective of the PLANCK mission is to provide a definitive characterisation of CMB temper- 
ature anisotropics, a good detection of E mode polarisation and possibly measure B mode polari- 
sation for the first time, see Tauber et al. [2009]. PLANCK is equipped with two instruments: the 
HFI is a bolometric array covering frequencies between 100 and 857 GHz Lamarre et al. [2009], 
while the LFI is a radiometric array covering the range of frequencies between 30 and 70 GHz, see 
Bersanelli et al. [2009]. 

Accurate measurements of the CMB fluctuations require a robust detection and removal of 
foreground emission. Data analysis strategies rely on their different frequency dependence in order 
to distinguish them from the CMB. Figure [I] shows the frequency dependence of dust, synchrotron, 
free-free diffuse emissions and the contribution from extragalactic sources, compared to the CMB. 
Thanks to its wide frequency coverage, PLANCK will be able to separate with high precision the 
different components by combining maps from different channels. 

Both radiometers and bolometers have a broadband response and their output is a convolution 
of the input spectra, weighted by the detector beam pattern, with their own frequency response. 
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Figure 1. The spectra of synchrotron, dust, and free-free foregrounds superimposed on the PLANCK chan- 
nels' bandwidths and the expected CMB signal at 1 degree scale (image from Efstathiou et al. [2005]) 

Good knowledge of the instrument's frequency response allows, if necessary, a fine correction 
to be applied at map level by an iterative process involving separation of the different components 
using their spectral indices, convolution of each foreground with the instrument frequency response 
and comparison of their sum with the observed map. 

Bandpass knowledge is also important for CMB polarisation measurements, the new frontier 
of precision cosmology that PLANCK is expected to open up. The HFI has polarisation capability at 
100, 143, 217 and 353 GHz thanks to Polarisation Sensitive BolometersRosset et al. [2009] (PSB) 
and all LFI channels are polarisation sensitive. Each LFI feed horn Villa et al. [2009] is connected 
to an orthomode transducer D'Arcangelo et al. [2009a] which splits the signal into two orthogo- 
nal linear polarisation components detected by two independent radiometers. Each radiometer is 
calibrated separately against the CMB dipole signal and polarisation measurements are performed 
combining the signal received by pairs of feed horns located symmetrically in the focal plane and 
following each other in the scan circle. For a complete description of the LFI focal plane see 
Bersanelli et al. [2009]. The main issue is that foregrounds have different spectral indices with re- 
spect to the CMB dipole which determined the calibration, therefore a spurious polarisation signal 
is generated when differencing channels with different bandpasses, producing leakage of the fore- 
grounds' total intensity into the polarised components. Using the radiometers' frequency response, 
it is possible to estimate the contamination of bandpass mismatch on CMB maps, as described in 
Leahy et al. [2009]. 

During qualification and performance tests, each active and passive unit composing LFI ra- 
diometers was characterised in terms of frequency response. For each unit we have measured 
return loss, insertion loss and (when applicable) gain and noise temperature as a function of fre- 



-2- 



quency over the entire bandpass. These hardware measurements are combined together by the LFI 
Advanced RF model, developed since 2001 Battaglia [2002], Franceschet [2004], which computes 
the bandpass of each channel of the integrated radiometer taking into account multiple reflections. 
For a complete review see Battaglia et al. [2009]). 

During the Flight Model (FM) test campaigns we attempted an end-to-end measurement of 
the radiometers' bandpasses on the integrated radiometer chains (see section |]), but this suffered 
from systematic effects in the measurement set-up. We have nevertheless been able to draw some 
qualitative conclusions from them. 

2. LFI radiometers' frequency response: basic equations and definitions 

The LFI gathers photons from the sky through 11 corrugated feed horns. Each horn is connected to 
an orthomode transducer (OMT) which splits the incoming signal into two orthogonally polarised 
components that propagate through two independent pseudo-correlation receivers that continuously 
compare the sky signal with a stable reference load at the temperature of ~4 K (Valenziano et al. 
[2009]). 

Each of the 22 receivers has cryogenic front-end module (FEM at ~20 K) and warm back- 
end module (BEM at ~300 K) radio frequency (RF) amplification stages (Davis et al. [2009], 
Artal et al. [2009], Varis et al. [2009]) connected by two composite waveguides D'Arcangelo et al. 
[2009b] and two square-law detectors that convert the RF power into a DC signal (figure 0). 




Figure 2. The 30 GHz Radiometer Chain Assembly (RCA) mounted on the support used during the cryo- 
genic test campaign on the individual receivers. 

In section ^7T] we describe the broadband response of the LFI radiometers. In section we 
shall briefly present the naming conventions used to identify the parts of the radiometers. 

2.1 Radiometer response 

Under the assumption of receiver linearity, the voltage output from each of the 44 output diodes 
can be represented mathematically by this equation: 

V ut =G-Pin J r V noise , (2. 1) 

where Pj n represents the input RF power, G is the photometric calibration constant and V no \ se is the 
noise contribution from the RF amplifiers and back-end electronics. 
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It is common to refer this noise contribution to the input, so that: 



Vout — G ■ (Pi n + Pnoise), (2.2) 

where P„ oise = V no ise/G. 

The power term Pnoise can be expressed in terms of noise temperature f no i se : 

j, 1 noise ,~ ~s 

* noise — ~, 7 ; y^-^J 

ksAv 

where kg is Boltzmann's constant, Av is the bandwidth, f no j se is the input load temperature that 
would produce an output signal V no j se in an ideal system, i.e. a system with no noise with V out = 

GPin- 

The noise temperature is directly related to the radiometer sensitivity by the radiometer equa- 
tion: 

a N = K Tin t^ -, (2.4) 

where: 

• is the Rayleigh- Jeans brightness temperature of the input source, 

• On is the root mean square of the white noise, 

• K is a factor which is dependent on receiver architecture and for the LFI receivers K = Vl 
Bersanelli et al. [2009], 

• t is the integration time. 

If we explicit the frequency dependence, Eq. flZTTl) reads: 

Voutiv) = G( V )-P in (v)+Vnoise(v), (2.5) 

where the function G(v), i.e. the effective gain versus frequency, is also referred to as the receiver 
bandpass. The relationship between integrated and frequency dependent quantities are described 
by the following equations: 



Vnoise = / V noise (v)dV, (2.6) 

Jo 

/•DO 

Pi„ = / Pm(v)dV, (2.7) 



G = (2.8) 

Pin 

During the PLANCK full-sky survey, the photometric calibration constant G will be routinely 
calculated from the modulation of the CMB dipole across the sky and then used to produce tem- 
perature mapsCappellini et al. [2003]. 

A priori knowledge of the bandshape G(v) is important for a detailed understanding of the 
data, particularly to ensure control of potential systematic effects that may be introduced by the 
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presence of foreground components both in the temperature and in the polarisation CMB maps. 
As evident from equation (|2.8|) , calibration is tightly bound to the spectral index, i.e. frequency 
dependence of the input source. CMB, free-free, dust and synchrotron foregrounds have different 
spectral indices and different couplings with each radiometric channel bandpass and in principle 
each would have a different photometric constant. For this reason, using the photometric constant 
derived from the CMB dipole introduces a second order systematic effect on these components. 
Knowledge of the instrument bandpasses can be used to reduce the effect by using an iterative 
process to separate different components Leahy et al. [2009]. 

The receiver bandpass, G(v), can be determined by performing two measurements per each 
frequency at different input power levels: 

2.2 LFI naming convention 

Each of the 1 1 assemblies connected to a feed horn is called a Radiometer Chain Assembly (RCA) 
and is identified with a number ranging from 18 to 28 according to the following scheme: 

• RCA numbers 18 through 23: 70 GHz, 

• RCA numbers 24 through 26: 44 GHz, 

• RCA numbers 27-28: 30 GHz. 

Each of the two radiometers connected to an RCA is labelled either with M is labeled either with 
M or S, whether it is connected to the main arm or the side arm of the orthomode transducer (OMT). 
Finally the four detector output diodes of each RCA are labelled and 1 if they are connected to 
M and 1 and 1 1 if they are connected to S. Therefore the label LF 1 2 5M- 1 indicates the detector 
1 in the radiometer M, i.e. connected to the OMT main arm of the 44 GHz RCA LFI25. 

For a detailed review of the LFI instrument refer to Bersanelli et al. [2009]. 



3. Bandpass estimation 

Our most accurate method to measure the LFI bandpasses is based on measurements of individual 
components integrated into the LFI Advanced RF Model (LARFM) to yield a synthesised radiome- 
ter bandpass. 

The LARFM is a software tool based on the Advanced Design System (ADS) circuit simulator 
by Agilent 1 . The measured frequency responses of the various subsystems (feed-OMT, FEM, 
BEM) are considered as lumped S-parameter components. 

Waveguides are simulated with an analytical model, in order to reproduce the effect of their 
temperature gradient 2 and the effect of standing waves caused by impedance mismatch at the in- 
terfaces between the FEM and BEM. This is because the 1.8-meter long waveguides were not 
measured at unit level in cryogenic conditions. The model provides accurate agreement with the 
measured waveguide response in the conditions of the test measurements (300 K). 

'Agilent Technologies, Inc http://www.home. agilent.com/agilent/product.jspx?nid=-34346. 0.00 
2 Recall that the LFI waveguides connect the 300 K BEM to the 20 K FEM 
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Figure 3. Model schematic for each LFI channel: feed horn and OMT are characterized by their S- 
parameters; FEM by gain and noise temperature; waveguides are simulated analytically; BEM by return 
loss, transfer function and integrated noise temperature. 



3.1 Model implementation 

The LARFM reproduces the frequency response of the 44 channels when looking at the input load. 
The overall layout of the model is sketched in figure 0. In this section we discuss the key points 
of each part that are relevant for bandpass response characterization. The detailed implementation 
can be found in Battaglia et al. [2009]. The frequency dependence of all parameters is implicitly 
assumed if not otherwise stated. 



3.1.1 Input load 

The input load is used to provide different input power in order to compute the spectral response. It 
is simply modelled as a resistor with tunable temperature; its output is an uncorrected noise signal 
independent of frequency whose power spectral density is: 

Pin = k B T in [W/Hz]. (3.1) 

3.1.2 Feed horn and orthomode transducer 

The feed horn and orthomode transducer are modeled as a single lumped two-port S-parameters 
component, characterized by scattering data taken at ambient temperature, see D'Arcangelo et al. 
[2009a]. Reflection scattering parameter (Sll and S22) measurements were performed on the sub- 
assembly comprising feed horn and OMT. Return losses were obtained by injecting a signal into 
the rectangular port facing the FEM. For the insertion scattering parameters (S21 and S12), the 
results of the OMT tests are representative of the whole assembly, feed horn losses being very low 
(typically < 1 dB). 

The 30 GHz OMT performance was measured between 26.5 and 40 GHz. In this case, all 
S parameters have been extrapolated down to 21 GHz, based on the 44 GHz OMTs which were 
measured over a very large band. The extrapolation was based on the normalization of each mea- 
surement at its own central frequency. 

It should be noted that our calculation of the radiometers' frequency response does not include 
the frequency dependence of the optical system, i.e., of the feed-telescope beam pattern. This has 
been discussed in detail by Sandri [2005], Sandri et al. [2009]. A full account of the instrument fre- 
quency response will require the combination of the optical and radiometric spectral dependences. 
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3.1.3 Front End Module 



The performance of the FEMs was measured during validation tests performed by the FEMs' pro- 
ducers: Jodrell Bank Observatory (JBO) for 30 and 44 GHz devicesDavis et al. [2009], and Ylinen 
for 70 GHzVaris et al. [2009]. Data are available both for FEM subcomponents, i.e. low noise am- 
plifiers, hybrids and phase shifters, and for the integrated FEMs. The performance of sub compo- 
nents is strongly influenced by connections and by mutual interaction, therefore it is more reliable 
to consider tests of the complete FEM. 

The FEMs' response is very linear and can be represented accurately by a linear model: 

PoutFEMiy) = G FEM (v) k B (7}„(v) + T noiseFEM (v)), (3.2) 

where: 

• Pomfem(v) is the output power per unit bandwidth; 

• Gfem(v) is the dimensionless gain factor; 

• TnoiseFEM is noise temperature, mainly due to the FEM HEMT amplifiers Davis et al. [2009]. 

Gain and noise temperatures were computed using the Y-factor method, which consists of 
fitting the linear model presented above with measurements at two different temperatures: 



„ , \ PoutThighiy) PoutTlowiy) ,~ 

Gfem{v) ~ MWv)-Tw(v)) ' (3 " 3) 

r(v) = p °«™* h{ y\ , 0.4) 

PoutTlow(V) 

T T high (v)-Y(v)T l0W (v) 

inoiseFEM\y) — Y{v) - 1 ' 

The front-end modules are modelled as amplifiers characterized by a spectral dependence of 
gain and noise temperature as measured in the FEM units. The measurements were obtained with 
the Y-factor method at the cryogenic temperature for each LFI radiometer FEM channel. The 
return loss parameters S 1 1 and S22 were inferred from the mean return loss measured at ambient 
temperature. 

All measurements covered the full radiometer bands. However, at 30 GHz, the out-of-band low 
frequency tail from 25 to 21.3 GHz, the WR28 waveguide cutoff, could not be measured directly. 
In this case we extrapolated the out-of-band profiles using a set of measurements performed on a 
larger frequency span of one of the FEMs with a Vector Network Analyser (VNA). 

3.1.4 Waveguides 

The 44 LFI composite waveguides are 1.5 to 1.8 meters long, depending on the channel, and 
they have a complex routing to allow integration of the HFI in the central portion of the focal 
plane. Furthermore, the waveguides support the thermal gradient between the ~ 20 K front end 
and the ~ 300 K back end, with three intercepts at the three V-grooves (see Tauber et al. [2009]). 
This complex thermo-mechanical layout led to the choice of measuring the waveguide responses 
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only at room temperature, and then the results were extrapolated to the results at nominal thermal 
conditions. This choice is justified by the fact that, while the cryogenic setup for a direct cryogenic 
measurement would be extremely complex and likely to be affected by systematic effects, the RF 
model of a waveguide is very well defined once its manufacturing characteristics are known. 

We implemented a software simulator using several ADS internal rectangular waveguide com- 
ponents connected in cascade with temperatures linearly interpolated from 20 to 300 K. 

This model was compared successfully to ambient temperature measurements and it is able to 
reproduce correctly the S-parameters of the waveguides. 

3.1.5 Back End Module 

The BEM provides amplification (about 25 dB) and detection of the RF signal, and outputs a DC 
voltage that is integrated in time and digitalized by the LFI Data Acquisition Electronics. Each 
of the four channels of a BEM contains an RF amplification stage based on HEMT Low Noise 
Amplifiers, a bandpass filter, a square law detector diode, and a DC amplifier. 

The overall frequency response of each BEM has been measured for each channel and is used 
as an input to the LARFM. In particular, each channel is implemented with an Sll S-parameter 
component (just modelling the return loss), a noise source and a function which computes the 
incoming power and multiplies it by the BEM gain. All the return loss S-parameters, the transfer 
function (RF input power to Volt output gain) were directly measured at each BEM. The noise 
source temperature was derived by integrating the measured noise temperature. 

A detailed description of the BEM devices and performance tests is available in Artal et al. 
[2009] for 30 and 44 GHz and Varis et al. [2009] for 70 GHz. 

3.2 Model output 

The LARFM computes the power due to thermal noise emitted by the input load at a given tem- 
perature in steps of 0. 1 GHz, then propagates it through all RCA components considering multiple 
reflections. As discussed in the previous section, feed horn/OMT, FEM and BEM responses are 
defined by experimental data at each frequency, while the waveguides response is simulated by the 
software model itself. 

To compute the bandpass, the simulation is run twice at different input load temperatures 
(2.7 K and 200 K) and equation Q2.9D is used. The concept of the procedure is similar to that of 
a swept source measurement, where a monochromatic power input is swept through the bandpass 
whilst recording the voltage output. 

As an independent check, we have earned out a swept source test of all RCAs in the cryogenic 
chamber. Before presenting the bandpass results (section |5|) we shall briefly present the end-to-end 
measurements in section @]. 

4. End-to-end cryogenic measurements 

The relative shape of the bandpass of each RCA was measured as part of the RCA test campaign 
by injecting a monochromatic signal into the input horn sweeping through the operational band and 
recording the DC output as a function of the input frequency, as shown in figure |]. 
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Figure 4. Schematic of the experimental setup for 30/44 GHz swept source tests: a monochromatic input 
microwave signal is injected into the input load. The input frequency changes with time across the radiometer 
band. The bandpass was computed from the recorded output voltage using Eq. 



The 30 and 44 GHz RCAs were assembled and tested in Thales Alenia Space in Vimodrone, 
while the 70 GHz RCAs were characterised at DA-Design Ylinen in Finland. The experimental 
setup was slightly different in the two cases. At 30 and 44 GHz, the input load was a cylindrical 
blackbody cavity 20 cm high and with a diameter of 20 cm, providing a return loss of about 60 
dB. The injector was an open-ended WR28/WR22 waveguide located in such a way that the signal 
was propagated into the RCA horn by the input load back wall. The results show, in general, good 
agreement with the bandpasses estimated by LARFM (see figures || and ||). 



RCA27 channel S-ll 



— LARFM 




Frequency[GHz] 

Figure 5. LFI27S-11 at 30 GHz. Comparison of LARFM bandpasses with end-to-end swept source mea- 
surements performed on WR28 nominal bandpass, between 26.5 and 40 GHz. 

However, two different kinds of systematic effect hampered the measurements: in many 44 
GHz channels a lack of power was observed in the high frequency portion of the spectrum; fur- 
thermore, in four cases, standing waves were produced as a consequence of poor alignment of the 
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RCA26 channel M-01 




Frequency[GHz] 

Figure 6. LFI26M-01 at 44 GHz. Comparison of LARFM bandpasses with end-to-end swept source mea- 
surements. 

polarised RF injection into the input load. 

At 70 GHz, the microwave generator was connected to a small input load placed directly in 
front of the horn, which yielded a return loss below the -30 dB requirement. This caused standing 
waves between the injecting horn and the RCA feed horn. 

In conclusion, while the end-to-end tests yielded a useful qualitative cross-check, these tests 
did not produce reliable results for a quantitative analysis. The LFI bands are determined by the 
LARFM combined approach already discussed. 

5. Results 

Figures |5] and |5] show all the LFI bandpasses obtained by the frequency response data of each 
radiometer unit assembled by the LARFM. The 70 GHz channels show a low bandpass ripple, of 
about 10 dB, which is within scientific requirements. The spike between 60 and 61 GHz, below the 
low frequency cut-off, is due to a systematic effect present in all the BEM gain measurements and 
caused by the test setup. We removed this range from the bandpasses made available at the Data 
Processing Center in order to avoid possible spurious effects and therefore the frequency coverage 
is 61-80 GHz. The high frequency cut-off is not well defined in most of the channels. 

The 30 and 44 GHz bandpasses show a more complex shape, driven by the BEM spectral 
response, but still within ±10 dB. The low frequency cut-off is always well defined, while the 
high frequency cut-off is not well defined in RCA 24 and 26. However, comparing with the high 
frequency cut-off of RCA 25, it is expected that the additional bandwidth is very low. Frequency 
coverage is 25-50 GHz for the 44 GHz channels and 21.3-40 GHz for the 30 GHz channels. 
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Figure 8. 30 and 44 GHz bandpasses, formatted as figure |3] 
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6. Conclusion 



We have presented the bandpasses of all the Planck-LFI radiometers as computed using the LFI 
Advanced RF Model based on the measured response of single radiometer units. Measurements of 
single components (feed-horns, orthomode transducers, front-end modules, back-end modules) are 
obtained with standard methods, typically employing a VNA test set, and provide highly reliable 
results, with precision of order 0.1-0.2 dB over the entire band. For the waveguides we have used 
simulated inputs that closely match the measured data. The composite bandpasses are estimated to 
have a precision of about 1.5 to 2 dB. 

We also attempted an end-to-end measurement of the RCA spectral response in the cryofacil- 
ity as an independent check. Unfortunately, these measurements suffered some subtle systematic 
effects in the test setup (standing waves at 70 GHz; polarisation mismatch and low frequency cov- 
erage at 30 and 44 GHz), preventing an accurate cross-check. However, the comparison shows a 
general agreement within limits of the test reliability and repeatability. 

In figure we show the normalised bandpasses for all the 44 LFI channels. The spectral 
responses exhibit a moderate level of ripple, within ±10 dB, dominated by the response of the back- 
end modules at 30 and 44 GHz, which reduces the radiometer effective bandwidth. The measured 
band shapes are in line with the scientific needs of the instrument and the knowledge achieved in 
these measurements allow us to analyse the impact of the radiometers' frequency response on the 
scientific outputs. 

The LFI Advanced RF Model is currently being ported to QUCS 3 , an open source microwave 
circuits simulator. This development is expected to lead to a more flexible and linear implementa- 
tion, see Franceschet [2004], and to facilitate the interface with external scripts to run simulations 
and export results in view of the forthcoming PLANCK data analysis. 
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